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ABSTRACT 

This paper presents the fourth extension to the Very Long BaseHne Array (VLBA) Calibrator Survey, con- 
taining 258 new sources not previously observed with very long baseline interferometry (VLBI). This survey, 
based on three 24 hour VLBA observing sessions, fills remaining areas on the sky above declination -40° 
where the calibrator density is less than one source within a 4° radius disk at any given direction. The share of 
these area was reduced from 4.6% to 1.9%. Source positions were derived from astrometric analysis of group 
delays determined at 2.3 and 8.6 GHz frequency bands using the Calc/Solve software package. The VCS4 
catalogue of source positions, plots of correlated flux density versus projected baseline length, contour plots 
and fits files of naturally weighted CLEAN images, as well as calibrated visibility function files are available 
on the Web at http : //gemini . gsf c . nasa . gov/ vcs4| 
Subject headings: astrometry, catalogues, surveys 



1. INTRODUCTION 

This work is a continuation of the project of surveying the 
sky for bright compact radio sources. These sources can be 
used as phase referencing calibrators for imaging of weak ob- 
jects with very long baseline interferometry (VLBI) and as 
targets for space navigation, monitoring the Earth's rotation, 
differential astrometry and s pace geodesy . The method of 
VLBI, first proposed by Matveenko et al.l lfT965) . allows us 
to determine positions of sources with nanoradian precision 
(1 nrad « 0.2 mas). Several catalogues were compiled com- 
bining observations under various programs. The catalogue of 
sources observed under geodetic programs from 1979 through 
2004, ICRF-Ext2 (Fe v et al. 2004). contains positions of 776 
sources. In addition to that, positions of 2247 sources were 
determined in the framework of the VLBA Calibrator Survey 
projec t: VCSl (Beaslev et al. 2002), VCS2 ( Fomalont et al. 
12003) and VCS3 (Petrov et al. 2005). Since 364 sources are 
listed in both the ICRF-Ext2 and the VCS catalogues, the to- 
tal number of sources for which positions were determined 
with VLBI is 2659. Among them, 2269 sources, or 85%, are 
considered as acceptable calibrators: they had at least 8 suc- 
cessful observations at both X and S bands and the semi-major 
axis of the eiTor ellipse of their coordinates is less than 25 nrad 
(5 mas). 

However, the sky coverage of these sources is not uniform. 
Successful phase referencing requires a calibrator within at 
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least 4° from a target source with a precise position and 
known source structure. The probability of finding a calibra- 
tor from the combined ICRF-Ext2, VCSl, VCS2 and VCS3 
catalogues within 4° of any target above -40° is 95.4%. In 
this paper we present an extension of the VCS catalogues, 
called the VCS4 catalogue, mainly concentrating on the other 
4.6% of the sky where the source density is the lowest and 
on the brightest sources with flat spectrum previously not ob- 
served with VLBI under geodesy and astrometry programs. 
Since the observations, calibrations, astrometric solutions and 
imaging are similar to that of VCS l-3 . most of the details are 
described bv lBeaslev et al] j200 2') and'P etrov eTaP (120051) . In 
section ^we assess an a priori probability of source detec- 
tion. In section ^we describe the strategy for selecting 412 
candidate sources observed in three 24 hour sessions with the 
very long baseline array (VLBA) based on analysis of prob- 
abilities of source detection which is a further development 
of the statist i cal ap proach for source selection introduced by 
iPetrov et al. l ^2005"). In section ^we briefly outline the ob- 
servations and data processing. We present the catalogue in 
|5] and summarize our results in |6] 

2. ASSESSMENT OF AN A PRIORI PROBABILITY OF 
SOURCE DETECTION 

Having unlimited resources one could try to observe all 
sources stronger than some limiting flux density. However, 
only sources with bright compact components can be detected 
with VLBI and may be useful for phase referencing or geode- 
tic applications. Kellermann, Paul inv-Toth. & Davisi (119681) 
first showed that the distribution of sources over spectral index 
a {F oc zy"*"") has two peaks: one near q = -I (steep spectrum) 
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Fig. 1. — The probability density distribution at 8.6 GHz as a function 
of spectral index. The thin lines show the distribution of the steep spectrum 
source population (left) and the flat spectrum the source population (right). 

and another near a = (flat spectrum). Later, it was confirmed 
that extended objects dominate in the steep spectrum popula- 
tion (see the review by Kellermann & Owen 1988). In com- 
pact regions, which have the synchrotron mechanism of emis- 
sion, the peak in the spectrum caused by synchrotron self- 
absorption has frequenc ies higher than 1 GHz due to their 
small size Thus, if the dominating emission 

comes from compact regions, the spectrum of the total flux 
density will be predominantly flat or inverted. 

Let's determine the probability of detecting with the VLBA 
a source with a given total flux density and given spectral 
index. The empirical probability distribution function of 
sources in spectral indexes can be approximated as 



P(a) = ( l-a)N(a,as, CTs) + aA^(Q!,Q;f, (Tf) , 



(1) 



where a is the fraction of the flat-source population to the 
total population, N(a,ao,(T„) is the normal distribution with 
maximum at a,, and dispersion a^. The first term Nsia) = 
N(a,as,crs) represents the probability distribution for steep 
spectrum sources, the second one N{(a) = N(a, af.af) — 
for flat spectrum sources. Using results of Mingalie v et alJ 
for the RATAN-600 simultaneous broad-band spec- 
tra survey at centimeter wavelengths of sources around the 
North Celestial Pole, as well as 2.7 & 5 GHz data from 
IWright & Otrupcek (1990) PKSCat90 catalog with Hmiting 
flux density at 2.7 GHz of 250 mJy, we have found the follow- 
ing estimates of distribution parameters: a = 0.15, (Tf = 0.30, 
Off = 0.15, (Ts = 0.28, as = 0.80. Figure^shows the distribu- 
tion. 

If a source has a spectral index in the range of [a,a + da], 
the ratio of the probability that it belongs to the flat spectrum 
population, Pf{a), to the probability that it belongs to the steep 
spectrum population, P^ia), is equal to the ratio of their prob- 
ability densities: 



Pi(a) aNf(a) 



Psia) il-a)Nsia) 



(2) 



Since we assume that a source either belongs to the flat 
spectrum population or to the steep spectrum population, 
Pf{a) + Ps(a) = 1. Then we immediately find 

aNfia) 



Piia) = 
Ps(a) = 



(I -a) Nsia) + aNfia) ' 

(l-fl)jVs(a) 
il -a) Nsia) + aNfia) 



(3) 



(4) 



The probability that a given source will have a ratio r 
of the correlated flux density /Vorr to the total flux den- 




FlG. 2. — The probability density distribution at 8.6 GHz of the ratio of the 
correlated flux density at the baseline projections equal to the Earth radius to 
the total flux density among the flat spectrum source population. 

sity Ffot is quite differ e nt for these two populations. Re- 
sults of 'Kovalev et al.' ("2005") of the 2 cm VLBA survey 
(Kellermann et al. 1998) and results of analysis of the VCS3 
observing campaign provided us estimates of /Von- at baseline 
projections equal to the Earth radius. The total flux of these 
sources was measured with the RATAN-600. The histogram 
over 500 sources after smoothing and normalization gives us 
an estimate of the empirical probability density distribution 
of the ratio of correlated flux density at baseline projections 
equal to the Earth radius to the total flux density among the 
flat spectrum source population (figure|2}. 

According to our previous experience in processing snap- 
shot VLBA Calibrator survey observations with integra- 
tion time 4 minutes and bit rate 128 Msamples/sec with 
good weather conditions, sources with correlated flux density 
greater than 60 mJy are reliably detected. We set the detection 
limit of the VLBA to f;„„ = 70 mJy in order to have a little al- 
lowance for possible pointing errors and higher than normal 
system temperature due to bad weather. The probability to 
have the correlated flux density greater than this limit, /v„„, 
i.e. the probability of detection for a source which belongs to 
the flat spectrum population and has total flux density F,o, is 

P'lim 



PiiFtot) = Cf 



where C/(r) is the cumulative probability 



function defined as 



Cf(r)=l - / Dfir)dr, 



(5) 



and where Dfir) is the probability distribution function of the 
ratio of the correlated flux density to the total flux density and 
r is the ratio of minimal detected correlated flux density to 
total flux density. The plot of Dfir) is shown in figure |2] and 
plot of Cf(r) is shown in figure|3l 

For steep spectrum sources the probability density function 
of the ratio of the correlated flux density to the total flux den- 
sity is not well known. According to our analysis of compact 
steep spectrum sources detected in VCSl-3 experiments, we 
approximate 



Dsir): 



•20, 
, 0, 



if r<0.05, 
if r>0.05. 



(6) 



Analogously, we introduce the function Cs(r) for the steep 
spectrum source population: 



Csir)=l - / Dsir)dr 



(7) 
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Fig. 3. — The cumulative probability function of the ratio of the correlated 
flux density to the total flux density among the flat spectrum source popula- 
tion: Cf{r,Fco,r/F,oi) = Pf(r < Fcorr/F,a,). 

Then, assuming the spectral index of a source is precisely 
known, the probability to be detected at 8.6 GHz for a source 
with spectral index a and total flux density F,, measured at 
frequency /„ is 

mJ„.a, = «a,C,(j-^||L_), 

Plots of the probability of being detected as a function of 
spectral index are presented in figure|3 

If no information about the spectral index is available, the 
estimate of the probability to detect a source is (see figure|5|l 

-oo 

P,(a) Q ( —, ]] da. (9) 

3. SELECTION OF CANDIDATE SOURCES FOR ZONES 
WITH LOW CALIBRATOR DENSITY 

As the first step of source selection we computed the sky 
source density on a 6' x 6' grid and identified areas with dec- 
lination > -40° with no calibrator within a 3° disk. 

As the second step we selected all sources from the NVSS 
catalogue with flux density greater than 100 mJy at L4 GHz, 
with galactic latitude \b\ > 1°.5 which either were not pre- 
viously observed with VLBI under astrometry and geodesy 
programs, or which were observed, but did not collect enough 
data in order to be counted as calibrators. We used the CATS 
database ( Verk hodanov et alJ 19971) to gather all available flux 
density measurements and then computed a spectral index 
for each source. Sources with spectral index less than -0.5 
were deselected. We also deselected sources with extrapo- 
lated total flux density at 8.6 GHz less than 100 mJy. Sources 
for which the spectral index could not be reliably determined 
were removed from the list of candidates if they had flux den- 
sities at 1 .4 GHz less than 200 mJy. We also removed pairs 
of sources if the distance between components was less than 
0°.083, since it usually indicates a source with extended struc- 
ture. Those sources, which according to NASA Extragalac- 
tic Database (NED), were listed as galactic objects, H II re- 
gions or planetary nebulae, were removed as well. Finally, we 



generated the input list of 2479 candidates. For 1216 objects 
we had estimates of their spectral indexes, and we computed 
a priori probability of their detection using expression (|8}. For 
other sources we computed this probability using expression 
It should be mentioned that many sources have variable 
flux density and spectral index. These variations were not 
taken into account in evaluation of the probability distribution 
and equation (|8}. This makes the estimates of the probability 
of detection somewhat less reliable. 

If all these sources could be observed, the mathematical ex- 
pectation of the area without a calibrator within a disk with 
radius of 4° would be 0.57%. For economical reasons we 
were forced to reduce the list of candidates. In order to find 
a subset of sources in such a manner that the area with the 
source calibrator density less than one source in a disk with 
radius 4° would be minimal, we adopted the following iter- 
ative procedure. Initially, we selected all sources. For each 
point on the grid with declination > -40° , which before the 
VCS4 campaign had no calibrators closer than 4°, we com- 
puted the probability of that that point will have no calibrator 
closer than 4° after VCS4 observations as 

Hij=l[(l-Pi(a,F))d(iJ,k), (10) 

k=\ 

where d(i,j,k) is 1, if the distance between the point on the 
grid and the ^h source is less than 4°, and otherwise. 
The sum 

R = Y,Hu (11) 

grid 

is the mathematical expectation of the total area which has no 
calibrator within 4°. 
Then at each step of iteration for each source we compute 

7?^ = ^ 7? t/(/, ^) / ( 1 - Pi(a,F)), i.e., the contribution of the 

grid 

kth source to R. The source with the minimal contribution is 
removed, the probability H/j is updated, and the procedure is 
repeated until source remained in the list. 

We generated the final list of candidates by concatenating 
two lists: 1) a list of 300 sources selected according to the 
iterative process described above; 2) a list of 100 sources 
north of declination -30° with the highest probability of de- 
tection — these are the strongest flat spectrum sources previ- 
ously not observed in geodetic VLBI mode. We also added 
12 sources which had been previously observed and detected 
under geodetic programs, but never imaged with the VLBA. 
Of them, 5 sources fell in the areas with no calibrators, and 
therefore were considered as belonging to list 1, the other 7 
sources were added to list 2. The mathematical expectation of 
the area which would remain with no calibrator after observ- 
ing the final list of 412 objects was 1.77%. 

The rational for including the list of the one hundred bright- 
est sources was that the area of the sky without a calibrator 
within a disk with radius 4° is reduced too slowly with an in- 
crease in the number of candidates beyond several hundred. 
For most geodetic applications having a bright calibrator a lit- 
tle bit further from the target is more important than having a 
calibrator closer, but much weaker. On the other hand, for ac- 
curate phase referencing, a relatively close weak calibrator is 
often preferable to a stronger calibrator which is much further 
away. So, some compromise is needed in order to improve the 
calibrator list. 
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Fig. 4. — Probability of being detected at 8.6 GHz as a function of spectral index for sources with flux density at 1.4 GHz 0.1 Jy, 0.2 Jy, 0.5 Jy and 1.5 Jy. 
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Fig. 5. — The probability of detection of a source with unknown spectral 
index as a function of total flux density at 1 .4 GHz. 



4. OBSERVATIONS AND DATA PROCESSING 

The VCS4 observations were carried out in three 24 hour 
observing sessions with the VLBA on 2005 May 12, 2005 
June 12, and 2005 June 30. Each target source was ob- 
served in two scans. The scan duration was 120 seconds for 
sources from the list of 107 bright objects, and 235 seconds 
for other target sources. The target sources were observed 
in a sequence designed to minimize loss of time from an- 
tenna slewing. In addition to these objects, 93 strong sources 
were taken from the GSFC astrometric and geodetic catalogue 
2 004f_astro^. Observations of 3-4- strong sources from 
this list were made every 1-1.5 hours during observing ses- 
sions. These observations were scheduled in such a way, that 
at each VLBA station at least one of these sources was ob- 

•^ Ihttp : //gemini . gsf c ■ nasa ■ gov/ solutions/ as tro| 



served at an elevation angle less than 20°, and at least one 
at an elevation angle greater than 50°. The purpose of these 
observations was to provide calibration for mismodeled atmo- 
spheric path delays and to tie the VCS4 source positions to the 
ICRF catalogue (Ma et al. 1998j). The list of troposphere cal- 
ibrators^ was selected among the sources which according to 
the 2cm survey results ( Kovalev et al. 2005) showed the great- 
est compactness index, i.e. the ratio of the flux density of an 
unresolved detail to the total flux density at VLBA baselines. 
In total, 505 sources were observed. The antennas were on- 
source 60% of the time. 

Similar to the pre vious VLBA C alibrator Survey ob- 
serving campaign (Petrov etal. 2005), we used the VLBA 
dual-frequency geodetic mode, observing simultaneously at 
2.3 GHz and 8.6 GHz. Each band was separated into four 
8 MHz channels (IFs) which spanned 140 MHz at 2.3 GHz 
and 490 MHz at 8.6 GHz, in order to provide precise mea- 
surements of group delays for astrometric processing. Since 
the a priori coordinates of candidates were expected to have 
errors of up to 30", the data were correlated with an accumu- 
lation period of 1 second in 64 frequency channels in order to 
provide extra- wide windows for fringe searching. 

Processing of the VLBA correlator output was done in three 
steps. In the first step the data were calibrated and fringed 
using the Astro nomical Image Processing System (AIPS) 
(lGreisenlll988h . In the second step data were imported to 
the Caltech DIFMAP package ( SheDherdl f!997h and images 
were produced using the optimized automated procedure orig- 
inally suggested by Greg Taylor ( Pearson et al. 1994). We 
were able to reach the VLBA image thermal noise level'* for 
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most of our CLEAN images jWrobel & Ulvestadll2005l) . Er- 
rors of our estimates of correlated flux density of sources 
brighter than 100 mJy are determined mainly by the accu- 
racy of amplitude calibratio n, which for the VLBA, accord- 
ing to Wrobel & Ulvestad ( 2005), is at the level of 5% at 
1-10 GHz. This estimate was confirmed by comparison of 
the correlated flux density with the single-dish flux density 
which we measured with RATAN-600 in March 2005 for 
slowly varying sources without extended structure. The meth- 
ods of single-dish observations and data processing can be 
found in Kovalev et al. ( 1999). In the third step, the data were 
imported to the Calc/Solve software, group delays ambigui- 
ties were resolved, outliers eliminated and coordinates of new 
sources were adjusted using ionosphere-free combinations of 
X band and S band group delay observables of the 3 VCS4 
sessions, 15 VCSl-3 experiments and 3976 twenty four hour 
International VLBI Service for astrometry and geodesy (IVS) 
experiments^ in a single least square solution. The boundary 
conditions which require zero net-rotation of new coordinates 
of the 212 sources listed as defined in the ICRF catalogue 
with the respect to their positions from that catalogue were 
imposed in order to select a unique solution of differential 
equations of photon propagation. 

In a separate solution, coordinates of the 93 well known tro- 
pospheric calibrators were estimated from the VCS4 observ- 
ing sessions only. Comparison of estimates of coordinates of 
these sources with coordinates derived from analysis of 3976 
twenty four hour IVS experiments provided us a measure of 
the accuracy of coordinates from the VCS4 observing cam- 
paign. The differences in coordinate estimates were used for 
computation of parameters a and b{&) of error inflation model 
in the form of \/ (a aY + b{5)^ , where tr is an uncertainty de- 
rived from the fringe amplitude signal to noise ratio using the 
error propagation law and 5 is declination. More detail s about 
analysis and imaging procedure can be found in Beasle Tet al.l 
( 12002) and Petrov et al. (2005). The histogram of source po- 
sition errors is presented in Figure |6l 



TABLE 1 

The number of observed and detected sources 
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Fig. 6. — Histogram of semi-major error ellipse of position errors. The 
last bin shows errors exceeding 9 mas. 
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Data 


# sources 


# detections 


# calibrators 


detection rate 


calibrators 










a priori a posteriori 




List 1 


305 


199 


117 


45 % 65 % 


38% 


List 2 


107 


93 


82 


84 % 87 % 


77% 


Total 


412 


292 


199 


57% 71% 


48% 



In total, 292 out of 412 sources were detected and yielded 
at least two good points for position determination. However, 
not all of these sources are suitable as phase referencing cal- 
ibrators or as targets for geodetic observations. Following 
Petrov et al. (2005) we consider a source suitable as a cali- 
brator if 1) the number of good X/S pairs of observations is 
8 or greater in order to rule out the possibility of a group de- 
lay ambiguity resolution error; and 2) position error before 
reweighting is less than 5 mas following the strategy adopted 
in processing VCS observations. Only 196 sources satisfy this 
calibrator criteria. It should be mentioned that our criterion 
for suitability as a phase calibrator is rather conservative, and 
sources which fail this criterion may still be useful for some 
applications. Among 292 detected sources, 34 were previ- 
ously observed but did not have enough observations before 
the VCS4 campaigns to be considered as calibrators. Of these 
34 objects, new VCS4 observations of 27 of them provided 
enough information to classify them as calibrators. 

Tabled shows statistics of the a priori and a posteriori de- 
tection rate. The a priori probability gives us a reasonable esti- 
mate of the detection rate. The sources from list 1 are weaker 
than the sources from list 2, and for 40% of them the inte- 
gration time of 4 minutes was insufficient to qualify them as 
calibrators. The a posteriori detection rate turned out greater 
than the a priori rate, because we set the detection limit to 
70 mJy when computing the a priori detection rate. In fact, 
we were able to detect sources with correlated flux density as 
low as 50 mJy with integration time 4 minutes. 

5. THE VCS4 CATALOGUE 

The VCS4 catalogue is listed in Table |2] The first column 
gives source class: "C" if the source can be used as a cali- 
brator, "-" if it cannot. The second and third columns give 
IVS source name (B1950 notation), and lAU name (J2000 
notation). The fourth and fifth columns give source coordi- 
nates at the J2000.0 epoch. Columns /6/ and 111 give inflated 
source position uncertainties in right ascension and declina- 
tion (without cos 8 factor) in mas, and column /8/ gives the 
correlation coefficient between the errors in right ascension 
and declination. The number of group delays used for posi- 
tion determination is listed in column 191. Columns /lO/ and 
/1 2/ give the estimate of the flux density integrated over entire 
map in Jansky at X and S band respectively. This estimate is 
computed as a sum of all CLEAN components if a CLEAN 
image was produced. If we did not have enough detections of 
visibility function to produce a reliable image, the integrated 
flux density is estimated as the median of the correlated flux 
density measured at projected spacings less than 25 and 7 MA 
for X and S bands respectively. The integrated flux density 
has a meaning of the source total flux density with spatial fre- 
quencies less than 4 MA at X band and 1 MA at S band filtered 
out, or by another words this is the flux density from all details 
of a source with size less than 50 mas at X band and 200 mas 
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TABLE 2 
The VCS4 catalogue 





Source name 


J2000.0 Coordinates 


Errors (mas) 




Correlated flux density (Jy) 






















8.6 GHz 


2.3 GHz 




3ass 


IVS 


lAU 


Right ascension 


Declination 


Aa 


A<5 


Corr 


#Obs 


Total 


Unres 


Total 


Unres 


Band 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 




0006-363 


J0008-3601 


00 08 33.661411 


-36 01 25.05213 


6.20 


32.64 


0.697 


8 


0.07 




0.22 


0.10 


X/S 


C 


0010+336 


J0012+3353 


00 12 47.382197 


+33 53 38.47157 


0.58 


0.84 


-0.508 


41 


0.16 


0.13 


0.07 


0.06 


X/S 


C 


0012+319 


JOO 15+32 16 


00 15 06.147414 


+32 16 13.30953 


0.33 


0.55 


-0.211 


86 


0.20 


0.11 


0.13 


0.09 


X/S 


c 


0021-084 


J0024-0811 


00 24 00.672734 


-08 11 10.04881 


1.08 


2.13 


-0.252 


31 


0.10 


0.08 


0.09 


0.07 


X/S 


c 


0024-114 


J0026-1112 


00 26 51.443027 


-11 12 52.42503 


0.96 


1.68 


0.219 


30 


0.11 


0.05 


0.19 


0.09 


X/S 


c 


0032+612 


J0035+6130 


00 35 25.310617 


+61 30 30.76144 


0.93 


0.60 


0.200 


71 


0.14 


0.06 


0.23 


0.07 


X/S 




0033-088 


J0035-0835 


00 35 46.250383 


-08 35 54.04258 


0.83 


1.57 


-0.240 


86 


0.13 


0.12 


0.06 


0.06 


X 


c 


0052-201 


J0054-1953 


00 54 32.948443 


-19 53 01.00202 


0.42 


0.83 


-0.091 


59 


0.16 


0.10 


0.20 


0.05 


X/S 


c 


0052-125 


J0055-1217 


00 55 11.782596 


-12 17 57.09709 


0.25 


0.46 


-0.220 


90 


0.24 


0.18 


0.21 


0.15 


X/S 




0057+101 


J0059+1022 


00 59 46.769108 


+10 22 40.42531 


16.93 


17.91 


0.797 


2 










s 



Note. — Tablel^is presented in its entirety in the electronic edition of the Astronomical Journal. A portion is shown here for guidance regarding its form and contents. Units of 
right ascension are hours, minutes and seconds, units of declination are degrees, minutes and seconds. 



at S band. Column /1 1/ and /1 3/ give the flux density of unre- 
solved components estimated as the median of correlated flux 
density values measured at projected spacings greater than 
170 MA for X band and greater than 45 MA for S band. For 
some sources no estimates of the integrated and/or unresolved 
flux density are presented, because either no data were col- 
lected at the baselines used in calculations, or these data were 
unreliable. Column /14/ gives the data type used for position 
estimation: X/S stands for ionosphere-free linear combina- 
tion of X and S wide-band group delays; X stands for X band 
only group delays; and S stands for S band only group de- 
lays. Some sources which yielded less than 8 pairs of X and 
S band group delay observables had 2 or more observations 
at X and/or S band observations. For these sources either X- 
band or S-band only estimates of coordinates are listed in the 
VCS4 catalogue, whichever uncertainty is less. 

In addition to this table, the html version of this catalogue 
is posted on the Web''. For each source it has 8 links: to a pair 
of postscript maps of the source at X and S-band; to a pair 
of plots of correlated flux density as a function of the length 
of the baseline projection to the source plane; to a pair of fits 
files with CLEAN components of naturally weighted source 
images; and to a pair of fits files with calibrated uv data. The 
coordinates and the plots are also accessible from the NRAO 
VLBA Calibration Search web-page^. 

Figure presents examples of naturally weighted contour 
CLEAN images as well as estimates of the correlated flux 
density versus projected spacings. The S band data for the 
source J1531-I-3430 is an example of one of the weakest 
sources which we still were able to image. It turned out to 
be an object with an inverted radio spectrum of the milliarc- 
second emission. The images of JOS 17-0520 show the com- 
pact structure with jet components which were detected on 
the level of several mJy/beam only. The source J0856-1105 
has the highest VLBI flux density at X band among all the 
new VCS4 objects and will be useful for geodetic applica- 
tions along with a few tens of other compact VCS4 objects 
with high flux density at VLBA baselines. 

6. SUMMARY 

The VCS4 Survey has added 258 new sources not pre- 
viously observed with VLBI. Among of them, 199 sources 
turned out to be suitable as phase referencing calibrators and 

* l^^^^A/ger^^^^^f^^i^^^oWvcs4|^^^^ 

' http : / /www, vlba ■ nrao ■ edu/astro/calib| 



as target sources for geodetic applications. The area with 
source density less than one calibrator within a disk of ra- 
dius 4° in any given direction on the sky with S > -40° was 
reduced from 4.6% to 1.9% (refer to figure |5] and table 
After processing the VCS4 experiments, the total number of 
calibrators has grown from 2268 to 2472. This pool of cali- 
brators was formed from analysis of 18 VLBI Calibrator Sur- 
vey and 3976 twenty four hour IVS astrometry and geodesy 
observing sessions. 

The strategy of source selection based on a priori probabil- 
ities of source detection developed in this study was success- 
ful. The conservative a priori estimate of detection rate was 
57%, while the a posteriori detection rate is 70 %. The a priori 
estimate of the area with the calibrator density less than one 
object within a disk of radius 4° was 1.8 %, while the a poste- 
riori value is 1.9 %. 

The sky calibrator density for different radii of a search cir- 
cle and for three different declination zones is presented in 
table |3l A further search for new calibrators to decrease the 
area with insufficient sky calibrator density would be consid- 
erably less efficient, because the majority of the flat spectrum 
sources in those areas have already been observed. The re- 
maining sources are the steep spectrum sources, or sources 
with unknown spectral index, or very weak objects. Accord- 
ing to figure |5] the lack of information about spectral index 
of potential candidates makes the detection rate of almost any 
sample about the same. 



TABLE 3 

The probability to find a calibrator at any given direction within a circle of a 
given radius in selected declination zones. 



Search circle 




Declination zone 




radius 


[-90°, -40°] 


[-40°, -30°] 


[-30°, +90°] 


4°.0 


63.6% 


92.2 % 


98.6 % 


3°.5 


55.4% 


86.0% 


96.0% 


3°.0 


44.1% 


76.4 % 


89.7 % 


2°.5 


33.2% 


56.6% 


77.9 % 


2°.0 


22.7 % 


46.6 % 


60.7 % 


1°.5 


13.6% 


29.4 % 


40.0 % 


1°.0 


6.3% 


14.2 % 


20.1% 


0°.5 


1.6% 


3.7% 


5.4% 



Note. — All sources from 3976 IVS astrometric and geodetic sessions 
and 18 VCSl, VCS2, VCS3 and VCS4 expeiiments with the VLBA which 
are classified as caUbrators are taken into account. 
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Fig. 7. — From top to bottom. Naturally weighted CLEAN image.s at S-band (2.3 GHz). The lowest contour levels (2 steps) on images are plotted at "clev" 
levels (Jy/beam), the peak brightness — "max" values (Jy/beam). The dashed contours indicate negative flux. The beam is shown in the bottom left corner of 
the images. Dependence of the correlated flux density at S-band versus projected spacings. The error bars were computed on the basis of correlated flux density 
scatter around mean value. Naturally weighted CLEAN images at X-band (8.6 GHz). Dependence of the correlated flux density at X-band versus projected 
spacings. 
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Fig. 8. — Source sky density before (left) and after (right) the VCS4 observing campaign. The area above -40° declination with calibrator source density less 
than one source within a 4° radius disk is shown in black. 



At the same time, some relatively bright flat spectrum 
sources remained in the zones where there is at least one cali- 
brator in a disk with radius 4°. A sample of 675 such sources 
was observed in July 2005 with the VLBA in a follow-up 
VCS5 experiment (Y.Y. Kovalev et al. in preparation). The 
addition of these sources to the VLBA Calibrator list will not 
affect the calibrator sky density with the search radius of 4°, 
but will increase the density of calibrators with smaller search 
circles, which will be beneficial for many applications, for in- 
stance, the VERA project (Hon ma et al. 200 3). 

New VCS4 detected sources were observed in 2005 with 
the Russian Academy of Sciences 600 meter ring radio tele- 
scope RATAN for measurement of their instantaneous 1-22 
GHz spectra (Y.Y. Kovalev et al. in preparation). These data 
will allow us to determine compactness of these sources and 
their accurate spectral indexes. Using this information one 
can further improve the estimate of the a priori probability of 
source detection with VLBI as a function of source flux and 
spectral index which is essential for source selection in future 
surveys. 
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